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The factorial moments of the pion multiplicity distribu- 
tions are calculated with HIJING and UrQMD and found to 
be independent of the pr range included, in contrast to recent 
simulations with the linear a model which leads to large en- 
hancements for pions with transverse kinetic energies below 
200 MeV. This supports the use of the ratio of the factorial 
moments of low and high pr pions as a signal of "new" physics 
at low momentum scales, such as the formation of disoriented 
chiral condensates. 



A major goal of relativistic heavy-ion collision experi- 
ments is to explore the phase diagram of hot and dense 
matter. In addition to the anticipated transition to the 
quark-gluon plasma phase in which the individual 
hadrons have dissolved into a chromodynamic plasma of 
quarks and gluons, it is expected that chiral symmetry 
will be approximately restored in the hot collision zone. 
Signals of this latter type of phase transition may arise 
from the subsequent non-equilibrium relaxation of the 
chiral order parameter which is expected to exhibit large- 
amplitude long-wavelength isospin-polarized oscillations 
around the normal vacuum configuration, often referred 
to as Disoriented Chiral Condensates ||. One expected 
consequence would be an anomalous broadening in the 
distribution of the neutral pion fraction. However, this 
observable poses serious practical challenges and in fact 
early DCC experiments carried out at CERN and Fermi- 
lab were not able to discern such a signal ||. Therefore, 
especially with the RHIC facility at Brookhaven National 
Laboratory now becoming operational, the need has in- 
tensified for theory to identify specific observables that 
may be particularly informative. In the present paper we 
address the pion multiplicity distribution. 

Recent simulations with the linear a model Q (see be- 
low) have suggested that the induced oscillations of the 
chiral order parameter amplifies individual pion modes in 
the soft part of the spectrum and, as a result, leads to en- 
hanced fluctuations in the multiplicity distribution of the 
produced pions. In that work, the multiplicity distribu- 
tions were analyzed in terms of their factorial moments, 



f m = (N(N - 1)(N - 2) • • • (N - m + 1)) 
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with N being the particle multiplicity in the specified ra- 
pidity and momentum region and m denoting the order 
of the moment (the average (•) is taken over a sample 
of events). We note that the first factorial moment is 
simply the mean multiplicity, /i = (N). A more conve- 
nient normalization is provided by the reduced factorial 
moments, 
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since these are all unity if the multiplicity distribution is 
of Poisson form. Thus, the deviation of the higher-order 
reduced factorial multiplicity moments provide a direct 
indication that non-poissonian fluctuations are present. 
In Ref. Q it was found that the reduced factorial mo- 
ments for soft pions (those with a transverse kinetic en- 
ergy below 200 MeV) were significantly in excess of unity, 
while those for the harder pions were remained consistent 
with unity. This feature suggests that such an analysis 
be made for the early RHIC data. 

However, while the occurrence of reduced factorial mo- 
ments in excess of unity is an indication of non-trivial 
processes, it is not necessarily a unique DCC signature. 
Thus, it is important to ascertain how specific this phe- 
nomenon is. For this purpose, we compare the reported 
results of the linear a model with those prediced by the 
two most commonly used event generators, namely HI- 
JING § and UrQMD [§. HIJING is expected to give a 
particularly good approximation to the particle produc- 
tion and fluctuations at high transverse momenta, while 
UrQMD includes a detailed treatment of resonance for- 
mation and decays. As both types of process lead to 
highly correlated pion production, one may expect these 
models to yield some enhancement of the multiplicity 
fluctuations. Therefore it is of interest to perform a quan- 
titative comparison of the results. 

Let us first briefly describe the calculations in Ref. [4|. 
Using a semi-classical treatment of the linear a model [7\, 
the chiral field <fr(r) = (<r(r), 7r(r)) was prepared to rep- 
resent a rod-shaped source with a bulk temperature of Tq 
(typically 240 MeV) and a radius of R Q (6-10 fm). The 
rod was then endowed with a longitudinal Bjorken scal- 
ing expansion and the classical field equation was solved 
numerically on a cartesian lattice, using the comoving 
coordinates (r, 77) in place of (t, z) (and starting from 
t = t = 1 fm/c), 
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Due to the imposed longitudinal expansion and, later on, 
the self-induced transverse expansion, the field ampli- 
tudes decrease rapidly. When sufficient decoupling has 
been achieved, a transverse Fourier resolution is made at 
each value of r\. From the resulting expansion coefficients 
of 7r(r) and -k(r) one may then extract the coefficients 
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which represent the probability amplitudes for finding 
a particle with rapidity 7/ and transverse momentum k. 
(Here fl± is the cross section of the box employed in 
the calculation, mk is the transverse pion mass, m'l = 
rn^ + fc 2 , and the last factor compensates for longitu- 
dinal scaling expansion.) Thus the expected number 
of pions within a certain rapidity interval is given by 

71^ = j dr)\xk ( 7 7)| 2 j f° r each transverse momentum k 
and each charge state j. (These quantities become con- 
stant after the decoupling has occurred.) The actual 

multiplicity was then selected from the associated 
Poisson distribution. The total multiplicity N (which is 
also Poisson distributed) can be obtained subsequently 
by adding up all the pions emitted within the specified 
phase space, 
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where fco denotes the maximum momentum of the soft 
pions. 

In order to obtain sufficient statistics, a number of in- 
dependent Bjorken rods were treated. Due to the ther- 
mal fluctuations, each such "event" has a unique initial 
field configuration so, consequently, the final states differ 
in detail. In particular, the expected multiplicities nj: 
fluctuate from one event to the next. Additional fluctu- 
ation arises from the subsequent sampling of the actual 
(integer) multiplicities n based on the values expected for 
a particular event n. While the latter statistical process 
is poissonian by design, the event-to-event fluctuation of 
n is generally not. Depending on the isospin orientation 
of the relaxing chiral order parameter, certain modes are 
amplified preferentially and this mechanism is the origin 
of the anomalous fluctuations reflected in the factorial 
moments, as discussed in Ref. 

In the present investigation, we consider central 
Au+Au events (impact parameter b < 3 fm) at the 
planned maximum beam energy at RHIC (100 GeV per 
nucleon) . In order to set the stage, we show in Fig. [I] the 
average multiplicity of positive pions with a transverse 
momentum below 200 MeV/c, as obtained with the HI- 
JING event generator ||. These soft pions constitute 
only a relatively small part of the total number of pions 
(the calculated mean px of midrapidity pions range from 
330 MeV (UrQMD) to 370 MeV (HIJING)). 
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FIG. 1. Rapidity distribution of soft positive pions 
(p T < 200 MeV), as generated by HIJING for a sample of 
central Au+Au events at RHIC (impact parameter b < 3 fm). 
The corresponding results for the negative and neutral pions 
are similar. 
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FIG. 2. The multiplicity distribution of positive pions in 
the mid-rapidity bin (\y\ < 0.5) for central Au+Au events at 
RHIC (b < 3 fm), as generated by HIJING (solid squares). 
The corresponding Poisson distribution having the same mean 
value is also depicted (open circles). 

We focus now on the multiplicity distribution in the 
mid-rapidity bin, \y\ < 0.5. Figure || shows the calcu- 
lated multiplicity distribution for positive pions of any 
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energy, as generated by HIJING for a sample of cen- 
tral Au+Au collisions. Also shown is the Poisson dis- 
tribution having the same mean multiplicity. The calcu- 
lated multiplicity distribution is significantly wider than 
the Poisson distribution. While this feature is due in 
part to the averaging over different impact geometries, it 
also reflects the presence of non-statistical components in 
the microscopic processes. The same feature still holds 
even if the HIJING analysis is restricted to the soft pions 
only. Moreover, similar features are also obtained with 
UrQMD for either grouping of the pions (not shown in 
the figure). Thus, it should be expected that the corre- 
sponding higher factorial moments will exceed unity. 

That this is indeed borne out is evident from Fig. [|, 
which shows the reduced factorial moments for both soft 
and hard pions as calculated with various models: HI- 
JING with and without jet quenching, UrQMD, and the 
results obtained in Ref. B for the Bjorken rods with the 
linear a model. It is evident that the latter results stand 
out: Whereas all calculations with HIJING and UrQMD 
yield a rather similar behavior, namely a gentle increase 
of F m with the order m, the linear a model leads to re- 
duced moments that remain close to unity for the hard 
pions while increasing rapidly for the soft pions. 
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FIG. 3. Reduced factorial moments F m for the multiplicity 
of positively charged pions emitted at mid rapidity in central 
Au+Au collisions at the full RHIC energy, as generated by 
HIJING with (squares) or without (circles) jet quenching and 
by UrQMD (diamonds). Also shown are the results obtained 
in Ref. [4] for idealized Bjorken rods (triangles). Solid symbols 
represent soft pions while open symbols represent hard pions. 

This qualitative difference in the behavior can be made 
more visible by considering the ratios between the re- 
duced moments for soft and hard pions, as shown in 



Fig. § While all the HIJING and UrQMD calculations 
predict a similar behavior for pions with low and high 
transverse momenta, the dynamical simulations with the 
linear a model yield a strong enhancement of the fluc- 
tuations in the number of \ow-px pions. In particular, 
although both HIJING (with or without jet quenching) 
and UrQMD produce some enhancement of the multi- 
plicity fluctuations (above a pure Poisson behavior), nei- 
ther one shows any distinction between soft and hard 
pions in this regard. Thus it appears that a compari- 
son of the factorial moments of the soft and hard pion 
multiplicity distributions may provide a useful observ- 
able which could indicate the presence of interesting dy- 
namics beyond what has been included in the standard 
event generators. In particular, a relative enhancement 
of the low-pT factorial moments may signal the forma- 
tion of disoriented chiral condensates and may therefore 
be used to gain experimental information on the global 
chiral properties. 
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FIG. 4. The ratios between corresponding values of the 
calculated reduced factorial moments F m for soft and hard 
pions, for the four cases displayed in Fig. ^. 

Let us finally emphasize some caveats associated with 
the present analysis. Although the two Monte-Carlo 
models used here successfully reproduce many features 
of particle production in hadronic and nuclear collisions, 
many issues regarding multiplicity fluctuations remain 
unresolved. In particular, the observed intermittency sig- 
nal associated with very small rapidity bins (Ay <C 1) 
H cannot be reproduced satisfactorily with the models 
employed here. However, for rapidity bins of the larger 
widths employed in the present investigation (Ay 1), 
the data of Ref. || can be reproduced within the ex- 
perimental error bars. As for the linear a model, it 
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is expected that the semi-classical treatment underesti- 
mates the enhancement (by a factor of two or more) || . 
Moreover, it applies only to baryon-free systems (hence 
our focus is on mid-rapidity pions) and current calcula- 
tional capabilities are limited to the relatively schematic 
Bjorken rod geometries treated in Ref. 0. Therefore, 
the results should be regarded as qualitative only. For- 
tunately, though, the conclusions presented here do not 
depend on the exact modeling but rely only on the gen- 
eral features of the dynamics. 

The present study was motivated by the recently re- 
ported finding that dynamical simulations with the lin- 
ear a model for idealized systems lead to strong enhance- 
ments in the multiplicity fluctuations for soft pions, while 
the harder ones display fluctuations of Poisson form. Our 
present investigation has shown that although both the 
HIJING and the UrQMD event generators also produce 
multiplicity fluctuations in excess of pure Poisson statis- 
tics, their magnitudes are relatively small (as compared 
to the DCC case) and, importantly, they do not depend 
on the pt range considered. These findings lend support 
to the adoption of this observable as an indicator of DCC 
formation. 

Thus, in conclusion, our analysis suggests that a strong 
relative enhancement of the multiplicity fluctuations of 
pions with low transverse momenta may provide a ro- 
bust DCC signal in high-energy nuclear collision experi- 
ments, such as those underway at RHIC. It is a special 
advantage that the basic observable, namely the number 
of pions in a given px range, should be readily obtainable 
experimentally. Thus, the suggested analysis qualifies as 
a Year-One task at RHIC. 
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